Numerical  Simulations  of  the 

Slingatron 

by  Gene  R.  Cooper,  Derek  A.  Tidman,  Mark  Bundy, 
and  Stephen  Wilkerson 

ARL-RP-30 


September  2001 


A  reprint  from  Proceedings  of  the  10th  U.S.  Army  Gun  Dynamics  Symposium, 

Austin,  TX,  23-26  April  2001. 


Approved  for  public  release;  distribution  is  unlimited. 


The  findings  in  this  report  are  not  to  be  constuied  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other 
authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  thereof 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 
to  the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5066 


ARL-RP-30  _ September  2001 


Numerical  Simulations  of  the 
Slingatron 


Gene  R.  Cooper,  Mark  Bundy,  and 
Stephen  Wilkerson 

Weapons  and  Materials  Research  Directorate,  ARL 

Derek  A.  Tidman 

ALCorp. 


A  reprint  from  Proceedings  of  the  10th  U.S.  Army  Gun  Dynamics  S3mposium,  Austin, 
TX,  23-26  April  2001. 


Approved  for  public  release;  distribution  is  unlimited. 


Abstract 


The  slingatron  mass  accelerator  is  described  for  several  track  configurations 
(shapes),  and  numerical  simulations  of  this  accelerating  mass  traversing  a  given 
track  configtiration  are  presented.  The  sled  is  modeled  as  a  point  mass  that 
interacts  with  the  slingatron  track  using  both  a  conventional  and  a  new  empirical 
velocity  dependent  friction  law.  The  closed  loop  circular  slingatron  was  found  to 
produce  high  maximum  sled  velocities  provided  flie  gyration  angular  speed  is 
always  increasing.  In  contrast  several  spiral  shaped  slingatron  tracks  reveal  that 
high  maximum  sled  velocities  are  obtainable  with  the  gyration  speed  held 
constant  In  fact,  a  shngatron  constructed  out  of  semi-circles  is  shown  capable  of 
generating  high  velocity  sleds  in  such  a  way  that  no  initial  sled  injection  is 
necessary.  Qioosing  the  proper  initial  gyration  phase  with  an  empirically 
determined  friction  model  allows  the  mass  sled  to  gain  ever-increasing  velocities 
when  placed  in  a  semi-circle  slingatron.  The  sled  bearing  pressure  and  its  total 
acceleration  are  examined  and  presented. 
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Abstract 

The  slingatron  mass  accelerator  is  described  for  several  track  configurations 
(shapes),  and  numerical  simulations  of  this  accelerating  mass  traversing  a  given  track 
configuration  are  presented.  The  sled  is  modeled  as  a  point  mass  that  interacts  with  the 
slingatron  track  using  both  a  conventional  and  a  new  empirical  velocity  dependent 
friction  law.  The  closed  loop  circular  slingatron  was  found  to  produce  high  maximum 
sled  velocities  provided  the  gyration  angular  speed  is  always  increasing.  In  contrast 
several  spiral  shaped  slingatron  tracks  reveal  that  high  maximum  sled  velocities  are 
obtainable  with  the  gyration  speed  held  constant.  In  fact,  a  slingatron  constructed  out  of 
semi-circles  is  shown  capable  of  generating  high  velocity  sleds  in  such  a  way  that  no 
initial  sled  injection  is  necessary.  Choosing  the  proper  initial  gyration  phase  with  an 
empirically  determined  fiiction  model  allows  the  mass  sled  to  gain  ever-increasing 
velocities  when  placed  in  a  semi-circle  slingatron.  The  sled  bearing  pressure  and  its  total 
acceleration  are  examined  and  presented. 
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Nomenclature 


A  = 

acceleration  vector 

A.= 

n  component  of  A 

Ai  = 

m  component  of  A 

D  = 

sled  diameter 

F= 

force  vector  acting  on  sled  F  =  |F| 

F//  = 

n  component  of  F 

Fx  = 

m  component  of  F 

f  = 

spin  frequency  of  gyration  vector 

f  = 

time  derivative  of  gyration  frequency  f 

i,j,k  = 

Cartesian  unit  triad 

L= 

sled  length 

M  = 

mass  of  sled 

m  = 

normal  unit  vector  on  concave  side  of  track 

n= 

normal  unit  vector  anti-parallel  to  the  track 

P  = 

bearing  pressmre 

R  = 

sled  radius  vector  |r|  =  R 

r  = 

gyration  vector  |r|  =  r 

s  = 

displacement 

t  = 

time 
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v= 

v= 

x  = 

y= 

a  = 

5  = 

e= 

r= 

x= 

p= 

<!>= 

w= 

{')= 

0= 

0= 
x  = 

Lxj= 


velocity  V  = 

relative  sled  velocity 
abscissa  of  sled 
ordinate  of  sled 

orientation  angle  =^-2;z'|_(^/2;7)J 

differences  between  upper  and  lower  semi-circle  radii 

lock-in  angle 

ratio  of  specific  heats 

radius  of  ciirvature 

coefficient  of  fiiction 

diameter  of  circle 

orientation  of  R 

orientation  of  r 

initial  value  of  y/ 

absolute  value 

d(  )m 

d(  )/dt 

average  value 

vector  product 
greatest  integer  <  x 
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Introduction 


A  mechanical  method  for  accelerating  a  mass  to  high  velocities  has  been  proposed 
by  D.  A.  Tidman^  using  a  device  called  the  slingatron^  Several  closed  loop  slingatron 
configurations  have  previously  been  examined.  For  these  configurations  the  accelerated 
mass  (called  a  sled  in  this  report)  interaction  with  the  slingatron  track  was  modeled  with 
magnetic  levitation^  or  with  a  mechanical  fiiction  force  that  is  proportional  to  the  normal 
force  exerted  by  the  track^  on  the  sled.  This  report  presents  numerical  simulations  of 
slingatrons  having  several  different  track  configurations.  The  sled  and  track  interaction 
continues  to  be  treated  as  a  normal  force  fnction  model  but  now  the  fiiction  coefficient, 
fi ,  is  either  a  constant  or  an  empirically  determined  fimction  depending  on  the  sleds 
velocity  relative  to  the  slingatron  track. 

These  simulations  reveal  that  the  mass  can  be  accelerated  to  very  high  velocities 
for  each  track  configuration  (track  shape)  examined  here.  Generally,  the  calculations 
show  that  spiral  shaped  slingatrons  will  produce  high  velocities  for  constant  gyration 
speeds  and  the  single  loop  circular  slingatron  will  produce  high  velocities  for  sufficiently 
high  gyration  acceleration.  We  emphasize  that  the  main  reason  this  happens  is  that  the 
sled  locks  into  a  constant  phase  angle  as  it  transverses  the  slingatron  track.  When  this 
occurs  the  sleds  velocity  is  optimally  increased  for  each  Itt  revolution  of  the  sled.  This 
will  continue  until  the  fiictional  force  dominates  the  component  of  the  Coriolis  force  that 
is  parallel  to  the  slingatron  track.  However,  there  are  cases  where  parameters  can  be 
chosen  so  the  parallel  component  of  the  Coriolis  force  is  always  greater  than  the  force 
due  to  fiiction.  In  such  cases  one  will  find  that  the  non-relativistic  sled  velocity  increases 
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indefinitely.  In  the  spiral  slingatrons  lock-in  occurs  almost  instantaneously  while  for  the 
single  loop  circular  slingatron  lock-in  occurs  after  a  sufficient  amount  of  time  has  passed. 

Theory 


Figure  1  Schematic  of  a  general  slingatron 


Figure  1  shows  a  sled  of  mass  M  moving  along  a  track  with  position  R  +  r  and 
corresponding  velocity  V  =  R  +  f .  The  sleds  velocity  is  assumed  large  enough  so  that  it 
stays  in  contact  with  the  track  at  all  times.  In  fact,  the  velocity  may  be  large  enough  to 
cause  mass  loss  due  to  abrasion  or  evaporation,  i.e..  M  9^  0 .  The  force  F  acting  on  the 
sled  can  be  written  as  F  =  Fj^m  +  F/^n  where  m  and  n  are  unit  vectors  pointing  normal 
(toward  the  concave  side  of  the  track)  and  anti-parallel  to  the  track.  Therefore, 
n  =  -R'/R'  and  m  =  n  x  k  for  unit  vector  k  pointing  out  along  the  normal  to  the  plane 
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of  the  track.  Calling  the  unit  vectors  along  the  abscissa  and  ordinate  i  and  j  allows  the 
momentum  equations,  see  Fig.  1,  to  be  written  as 


(Mx  +  Mx)=  -(sin((t))(F^R'-  F,r)+  cos((|)Xf^R  +  F,,R'))/VR'^  +  R^ 

(My  +  My)=  -(sin((j)XFj.R  +  F//R')+  cos((t)XF//R  -  F^R'))/VR'^  +  R^  ^  ^ 

with  Cartesian  coordinates  (x,y)  in  the  (i,  j)  plane  which  are  constrained  so  that 

X  =  R  cos(^)  +  r  cos(^^a) 
y  =  R  sin(^) + r  sin(y/) 

The  parallel  force  ¥„  is  assumed  to  be  caused  by  friction  and  therefore  the  standard 
friction  model  is  represented  as 

F„  =  s  V  =  |^|V5JT7  (3) 

for  which  V  is  the  sled  velocity  relative  to  the  track.  Now  solving  Eqs.  (1)  for  Fj^  and 
F//  and  then  including  Eq.  (2)  and  Eq.  (3)  results  in  the  following  differential  equation 


V_  <^^(//(R^-RR'^  +  2R^")+RtR  +  R”)) 

R"+R'^ 

^  ((M(^y  +  //y}-^)+M;y)sin(y/--^)+(M^y^-//(My^  +  My;-))cos(^y-^))rR' 


M(R'+R'2) 

{{uiy/^  -  juMy/)sm(y/  +  +M^y  )  cos(y^-^))rR' 

M(R'+R'2) 


M 


Equation  (4)  is  general  enough  to  account  for  many  slingatron  configurations  as  long  as 
the  motion  takes  place  in  the  i,  j  plane.  The  configurations  are  obtained  by  specifying 

the  <j>  dependence  of  the  vector  R(^)  giving  the  shape  of  the  slingatron  track.  The 
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jfriction  coefficient, //(v),  the  time  dependent  mass,  M(t),  and  the  gyration  angle,  ^^/■(t), 
must  also  be  given  in  order  to  nmnerically  integrate  Eq.  (4) 

Circular  Slingatron 

This  report  considers  several  possible  track  configurations  all  of  which  treat  the 
length  of  the  gyration  arm,  |r|  =  r ,  as  a  constant.  The  first  consideration  is  the  circular 
slingatron'’^  defined  by 

|R(t»)|  =  constant .  (5) 


In  general  we  have  foimd  that  any  slingatron  will  optimally  accelerate  a  sled  to  high 
velocities  whenever  the  lock-in  angle,  6  =  \i/  -(p  ,\s  very  close  to  constant  i.e. 


0  w  0  hence 

Iff 


(6) 


Using  the  constraints  given  by  Eq.  (5)  and  Eqs.  (6)  simplifies  Eq.  (4)  to  the  following 
expression 


^ 

^  _  r(sin(^)  -  jj.  cos(^))  -  /IR  (7) 

r(^sin(0)+cos  ie)V^ 

for  constant  n  and  M.  This  is  easily  integrated  to  give  us  the  sleds  velocity 
y  =  as 


l^( 

'4 

/2rRcos(0)+  +  r^ 

l-<i(0)bt 

(8) 
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A  plot  of  Eq.  (8)  is  shown  in  Fig.  2  for  0  =  njA  ^^(o)  =  0  and  4>(o)  =k.  One  can  see  that 

the  non-relativistic  velocity  V  becomes  infinite  at  the  time  t^  =l/^^(o)b  whenever  we 
force  6  to  be  constant^ 


Fig.  2  Velocity  vs.  time  for  the  circular  slingatron 

The  values  of  t„  change  rapidly  in  regions  where  ^  « ;rn,  n  =  0,1,2"-  and  therefore  a 
plot  of  t„(^)  is  given  in  Fig.  3. 


Circular  Slingatron  r/R=0.2  f  (0)=it  ^=0.005 


Fig.3t.(e)  VS.  time  for  the  circular  slingatron 

Equation  (7)  shows  that  b  =  0  whenever 
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^  =  ±2  tan  '  {j{a^  +aj{a- 1)//)±  ;r  n,  n  =  0,1,2 

(\  ^  ' 
—  \)jj^  +  a  >0 

These  results  used  in  Eq.  (8)  produce  two  time  independent  velocities  given  by 

v  =  |iz^(o)|VR'-i]?  +  a^  +  ijy/ +1 .  A  small  amount  of  algebra  reveals  that 

constant  velocities  occur  when  the  normal  force  becomes  so  large  that  friction  prevents 
the  sled  from  accelerating  along  the  circular  track.  Hence,  the  sleds  velocity  stays 
constant. 

Instead  of  demanding  the  constraints  given  by  Eqs  (6)  lets  now  consider  a  circular 
slingatron  with  an  accelerating  angular  gyration  speed  such  that 

y^  =  y^o+2^ft  +  7r{t^.  (10) 

Numerically  integrating  Eq.  (4),  subject  to  Eq.  (10),  produces  the  velocity,  with  initial 
conditions  ^(O)  =  0  and  ^{o)  =  2ni ,  plotted  in  Fig.  4.  Here  we  see  that  high  velocities 
are  obtained  for  f  5^  0 . 


Fig.  4  Velocity  vs.  time  for  the  driven  circular  slingatron 
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The  corresponding  lock-in  angle,  0  =  ij/-^ ,  found  in  Fig.  5  becomes  relatively  constant 
when  sled  velocities  become  large. 

R=25m  r=1m  f=50cps  f'=10cps^  %=7t/4  p=0.005 

8.0E-01 

.e-  4.0E-01 

>  O.OE+00 
II 

®  •4.0E-01 
-8.0E-01 

0.0  1.0  2.0  3.0  4.0 

Time  (s) 

Fig.  5  Lock-in  angle  vs.  time  for  the  driven  circular  slingatron 
Even  though  high  sled  velocities  can  be  obtained  with  a  circular  slingatron  it  is 
encumbered  with  the  difficulty  of  creating  an  easy  exit  port  for  the  high-speed  sled.  To 
alleviate  this  mechanical  difficulty  we  will  next  examine  possible  open  loop  spiral 
slingatrons. 

Archimedes  Spiral  Slingatron 

We  now  focus  our  study  on  configurations  that  are  not  closed  loop  slingatrons.  In 
particular,  we  examine  spiral  shaped  tracks  having  open  ends,  which  have  the  advantage 
of  an  easy  exit  for  an  accelerated  sled  entering  fi*ee  flight.  Following  Tidman^  we  first 
examine  an  Archimedes  spiral  shaped  slingatron  track  given  by 

R  =  Ri+r(^sin({«/'o)  (11) 

for  which  the  constant  Rj  is  the  initial  radius  and  spacing  between  adjacent  spirals  is 
determined  by  sin(^Q).  The  gyration  speed  is  now  held  constant  so  that  the  expression 
given  in  Eq.  (10)  is  replaced  with 
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y/  —  y/Q+lTuft . 


(12) 


Tidraan^  presents  a  first  order  analysis  for  a  sled  with  M  =  0  and  jj.  held  constant. 
Putting  Eq.  (11)  into  Eq.  (4)  and  keeping  first  order  terms,  i.e.  terms  first  order  in 

,  leads  to  the  simple  equation  with  O^y/ 

V7V  =  (rsin(^))/R(^^J)-//.  (13) 

Integrating  this  result  shows  that  the  average  increase  in  velocity,  AV ,  for  the  sled  during 
a  In  revolution  of  a  single  spiral  leads  to  the  following  expression  for  averaged 
quantities 

AV/cycle  «  2^{Y)  [(r  sin(0))/(R(^))  -  //].  (14) 

With  this  motivation  we  now  place  Eqs.  (10-11)  into  Eq.  (4)  and  follow  with  numerical 
integration  to  produce  the  velocity  plot,  V(t)  =  ^Jx^  +y^ ,  found  in  Fig.  5.  The  initial 

conditions  are  ^(o)  =  0 ,  <^(o)  =  2;r  f  and  the  sled  mass  plus  the  coefficient  p  are  both 
held  constant. 


Fig.  6  Velocity  vs.  time  for  the  Archimedes  spiral  slingatron 


We  immediately  see  that  the  maximum  velocity  has  the  same  order  of  magnitude  as  that 
found  for  the  circular  slingatron,  see  Fig.  2,  but  does  so  in  far  less  time  and  requires  no 
acceleration  of  the  gyration  arm.  The  radius  changes  so  that  the  lock-in  angle,  d , 

stays  relatively  constant  until  the  sled  has  reached  its  maximum  velocity  as  shown  in 
Fig7. 


Fig.  7  Lock-in  angle  vs.  time  for  the  Archimedes  spiral  slingatron 
Inspecting  Eq.  (14)  one  can  infer  that  the  maximum  velocity  occurs  where  the  frictional 
force,  between  the  sled  and  track,  balances  the  parallel  component  of  the  Coriolis  force  as 
the  sled  moves  along  the  spiral  track.  A  plot  presented  in  the  following  section  will 
provide  numerical  confirmation  of  this  conclusion. 

Semi-Circle  Slingatron 

The  next  slingatron  configuration  examined,  is  constructed  using  a  sequence  of 
semi-circles^  having  increasing  radii,  see  Fig.  8. 
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Fig.  8  Schematic  of  the  semi-circle  slingatron 

We  note  this  slingatron  continues  to  have  the  sled  exit  advantage  found  for  the 
Archimedes  slingatron.  However,  this  track  has  the  added  advantage  in  that  large  sled 
velocities  can  be  obtained  without  the  need  of  initial  sled  injection,  i.e.  ^(o)  =  0  and 

^(O)  =  0 ,  provided  the  slingatron  has  been  scaled  up  to  a  sufficiently  large  size.  All  of 
the  previous  slingatrons  plus  the  small-scale  high  frequency  version  of  the  present 
slingatron  require  initial  sled  injection,  ^(0)9^0,  in  order  to  obtain  large  maximum 
velocities.  Because  of  these  advantages  a  more  thorough  investigation  including  sled 
pressure  loads  and  an  empirical  friction  model,  incorporating  mass  loss  M<0,  is 
presented.  According  to  Fig.  8  the  semi-circles  in  the  upper  half  plane  are  described  by 


R  =  Rj  +2  5L<|)/27rJ 
R'  =  0 
R'  =  0 

6  =  7rrsm(vj/o) 


(15) 


and  applying  the  law  of  cosines  to  this  geometry  we  find  for  the  semi-circles  in  the  lower 
half  plane  having  center  coordinates  at  (<5‘,0)  are  described  as 


13 


(16) 


R  =  <5 cos(^)+  ^(Rj  +2  5\<l>l27i:\+5j 
p,_  <?Rsin(^) 

(Jcos(^)-R 

p,_R^^+2«?  R^sin((»)+JR  cos(^) 

^cos(^)-R 

Putting  Eq.  (15)  and  Eq.  (16)  into  Eq.  (4)  followed  with  numerical  integration  leads  to 
sled  velocities  shown  in  Fig.  9. 


Fig.  9  Velocity  vs.  time  for  the  semi-circle  slingatron 

This  high  frequency  small-scale  slingatron,  i.e.  f  =  50cps  and  Rj  =  100cm,  is  subjected 

to  the  conditions  <f>{0),  <^(o)  =  2;r  f ,  M  =  0 ,  and  //  =  0.005 .  Figure  10  tells  us  again  that 

the  lock-in  angle  6  diverges,  from  nearly  a  constant  value,  after  the  time,  t ,  where  the 
sled  reaches  its  maximum  velocity. 
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Fig.  10  Lock-in  angle  for  the  semi-circle  slingatron 

An  important  simplification,  when  constructing  a  semi-circle  slingatron,  occurs 
when  we  consider  a  large-scale  version,  i.e.  f  =  llcps  and  Ri=400cm,  of  this 
slingatron.  For  these  cases  we  can  generate  large  maximum  velocities  without  the  need 
for  an  initial  sled  injection  velocity.  An  example  of  this  is  depicted  in  Fig.  1 1  such  that 

^(o)  =  0 ,  <^(o)  =  0  and  the  sled  mass  as  well  as  the  friction  coefficient  are  still  held 
constant. 


Fig.  11  Velocity  vs.  radius  for  the  large-scale  semi-circle  slingatron 

This  plot  has  the  time  axis  replaced  by  radial  distance  to  indicate  the  size  of  this  large- 
scale  example.  In  contrast  we  note  that  the  small-scale  version,  of  this  slingatron,  still 
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requires  a  non-zero  initial'  sled  velocity,  ^(o)  5^  0 ,  in  order  to  obtain  a  significantly  larger 
maximum  sled  velocity. 

To  address  a  more  realistic  slingatron  one  should  consider  the  interaction  between 
the  sled  and  the  slingatron  track.  Therefore,  we  will  assume  that  the  sled  is  covered  with 
material  that  abrades,  M  <  0,  as  it  moves  along  the  track,  in  such  away  that  the  removed 
mass  acts  as  a  lubricant.  As  the  velocity  increases  the  abraded  particles  may  become  a 
liquid  bearing  and  at  still  greater  velocities  the  liquid  will  eventually  evolve  into  a  gas  or 
even  a  plasma.  To  date  only  a  preliminary  experimental  investigation  of  the  friction 
coefficient,  ;/(v),  dependence  on  the  sleds  relative  speed  V,  has  been  carried  out  by 
Tidman^.  The  preliminary  data  for  a  lexan  sled  indicate  that 


V 


(v)  = 


0.152 
ll  + 3.16V 


=  ^VR''+R^ 


(17) 


where  V  is  measured  in  Km/s,  ranging  up  to  a  maximum  velocity  of  2.0  Km/s.  The 

corresponding  value  of  M  is  still  very  suspect  but  the  preliminary  data  gives  the 
following  estimation 


Ti/r 

50;dl(<z))‘ 

Tidman  has  suggested  that  if  the  velocity,  and  therefore  the  bearing  pressure  between  the 
track  and  the  sled,  is  large  enough  to  cause  the  abrading  mass  to  gasify  one  might  then 
assiime 
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(19) 


v,  jlMV 

Kr-i)RW 

forlexan/  » 1.25 

The  two  estimates  found  in  Eqs.  (18-19)  are  plotted  in  Fig.  12  while  assuming  the  friction 
coefficient  given  by  Eq.  (17)  remains  true  for  all  velocities. 


Fig.  12  Sled  mass  vs.  time  for  the  semi-circle  slingatron 
Figure  13  shows  the  velocity  results  for  a  constant  mass  sled  plus  the  two  cases  where 
M  0  given  by  Eqs.  (18-19).  Since  there  is  no  discemable  difference  between  the  three 
cases  we  will  select  M  =  0  for  the  remainder  of  this  paper. 


Fig.  13  Velocity  comparision  vs.  time  as  functions  of  M  and  M 
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This  selection  should  not  be  taken  as  a  general  conclusion  since  the  present  results  stem 
from  two  simple  M  models.  All  that  we  can  claim  here  is  that  the  dynamics  represented 
by  Eq.  (4)  are  weekly  dependent  on  our  two  M  models.  A  more  in-depth  study  of  the 
bearing  physics  may  reveal  that  M  cannot  be  ignored  in  a  more  elaborate  friction  model 
and  therefore  M  must  be  included  in  the  momentum  equations. 

To  demonstrate  the  importance  of  the  friction  model  we  next  compare  our 
empirical  model,  Eq.  (17),  to  the  constant  friction  model  n  =  0.005 .  Figure  14  has  an 

example  comparison  of  the  sled  velocity  using  initial  conditions  ^(o)  =  0  and  <^(o)  =  0 
for  these  two  models. 


Fig.  14  Velocity  comparison  vs.  time  as  a  function  of  p 
We  see  that  Eq.  (17)  and  the  initial  phase  xj/^  ^Tcje  can  lead  to  classically  unlimited 
velocities  imparted  to  the  sled.  This  is  an  example  force  due  to  friction  never  becomes 
large  enough  to  overcome  the  parallel  component  of  the  Coriolis  acceleration. 

To  gain  insight  into  the  magmtudes  of  the  bearing  pressure  exerted  by  the  sled  on 
a  semi-circle  slingatron  one  can  examine  a  typical  right  cylinder  shaped  lexan-coated  sled 
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having  length  L  and  diameter  D .  Assiiming  that  the  ablating  material  exerts  a  uniform 
pressure  along  the  half  lateral  surface  of  the  sled  facing  the  slingatron  track,  one  will  find 
the  pressure  P  is  easily  calculated  to  be  P  =  Fj^/DL .  An  example  of  bearing  pressure, 

for  a  sled  with  initial  mass  M  =  1 .0  x  1 0^  Kg ,  plotted  as  a  function  of  time  is  presented  in 
Fig.  15. 


Fig.  15  Sled  bearing  pressure  vs.  time  for  semi-circle  slingatron 
Calculating  the  sleds  acceleration  along  the  unit  vectors  m  and  -  n  while  the 
traversing  the  slingatron  track  may  also  prove  useful  to  the  designer.  Remembering  that 
the  normal  component  of  the  acceleration  is  always  pointing  toward  the  concave  side  of 
the  track  we  can  find  its  magnitude  from  Aj_  =  .  An  example,  for  the  large-scale 

shngatron,  is  shown  in  Fig.  16. 
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Fig.  16  Normal  component  of  acceleration  vs.  time  for  the  semi-circle  slingatron 

In  a  similar  fashion  the  parallel  component  of  acceleration.  A,,  =  -s,  is  plotted  in 

Fig.  17  where  for  comparison  we  also  plot  |Fj^|/M  . 


Semi-Circle  Slingatron 
Ri=4mr=2mfi=11cps'Fo=^t/4  p=p(v) 
D=76cm  L=434cm 


Fig.  17  Comparison  of  parallel  component  of  acceleration  and  friction 
force/weight  vs.  time 

Comparing  Fig.  16  with  Fig.  17  verifies  that  the  time  t  where  the  maximum  velocity  is 
found  is  also  the  time  where  fiiction  force  per  unit  mass  equals  A„ ,  thus  substantiating 
the  conclusion  mentioned  above. 

Examining  the  curvature  of  the  semi-circle  slingatron  shows  that  R"  is  not 
continuous  at  the  points,  (f)  =  n;r,  n  =  0,1,2 •••,  where  the  upper  half  plane  semi-circles 
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meet  the  lower  half  plane  semi-circles,  see  Eqn.  (15-16).  However,  these  discontinuities 
of  the  sleds  acceleration  are  not  severe  enough  to  cause  numerical  integration  problems 
of  Eq.  (4). 

Further  Illustrations 

Before  leaving  the  topic  of  spiral  shaped  slingatrons  we  will  give  results  for  two 
other  familiar  spiral  shaped  tracks.  Both  of  these  configurations  have  the  non-zero  initial 

velocity,  ^(o)  =  0  and  ^(o)  =  2;r  f .  The  first  is  the  parabolic  spiral  given  by 

R(^)  =  Rj  +  r  sin(^«/o  )  (20) 

A  velocity  plot  for  this  is  given  in  Fig.  18 


Fig.  18  Velocity  vs.  time  for  the  parabolic  spiral  slingatron 

for  which  the  fiiction  coefficient  //  is  modeled  using  Eq.  (17).  The  last  case  considered 
is  the  logarithmic  spiral  given  by 


R(^)  =  Rj  +  r  sin(^i^o)ln(^) .  (21) 

The  velocity  for  this  spiral  is  presented  in  Fig.  19  where  again  the  coefficient  n  is 
calculated  using  Eq.  (17). 
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Ri=400m  r=200m  f=11cps  'Po=7:/4  m^h(V) 

4.E+03 
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>  O.E+00 

0  5  10  15  20 

Time  (s) 

Fig.  19  Velocity  vs.  time  for  the  parabolic  spiral  slingatron 

One  can  see  from  the  last  two  plots  that  a  wide  variety  of  choices  for,  R(^^),r,vJ'  andf , 

are  possible  to  obtain  large  maximum  sled  velocities.  These  choices  offer  considerable 
flexibility  when  faced  with  design  constraints  for  a  spiral  shaped  slingatron. 

Conclusion 

The  closed  loop  circular  slingatrons  require  the  gyration  arm,  r,  to  accelerate, 
y/  >  a,  in  order  for  the  mass  sled  to  reach  substantial  maximum  velocities.  As  time 
progresses  the  sled  is  optimally  accelerated  as  indicated  by  the  lock-in  angles  6 
approaching  nearly  a  constant  value.  Experience  from  generating  computer  simulations, 

for  the  closed  circular  slingatron,  reveal  that  initial  sled  velocities  close  to  ^(o)  =  2;rf 
are  necessary  in  order  to  gain  substantial  sled  acceleration. 

The  spiral  slingatrons  are  able  to  produce  sufficiently  large  maximum  sled 
velocities  using  only  constant  angular  gyration  speeds,  \j/  =  0.  All  of  the  small-scale 

spiral  slingatrons  require  initial  injection  velocities  in  the  neighborhood  of  (/i(o)  =  2;r  f 
in  order  to  gain  large  maximum  speeds.  However,  the  large-scale  version  of  the  semi¬ 
circle  slingatron  has  an  additional  favorable  feature  in  that  it  can  generate  large  maximum 


1 

Logarithmic 

Spiral 

Slingatron 

t-:  ,  , _ 1 
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velocities  with  zero  initial  injection  sled  velocities.  This  will  make  such  slingatrons 
mechanically  easier  to  build.  Examination  of  the  bearing  pressure  and  the  corresponding 
magnitude  of  the  sled  acceleration  give  insight  into  the  stress  levels  that  the  sled  and 
spiral  track  must  endure. 

The  empirical  friction  model  shows  that  the  maximum  velocity  is  very  sensitive  to 
Fuorthermore,  we  found  unlimited  sled  velocities  for  a  proper  choice  of  the 
parameter  y/^ .  There  are  a  variety  of  choices  regarding  the  spiral  shape,  R(^i),  as  well  as 
the  parameters  andf  that  one  can  make  in  order  to  produce  a  large  range  of 

maximum  sled  velocities.  Even  though  the  results  given  here  strongly  suggest  that  spiral 
shaped  tracks  are  the  most  desirable,  the  optimal  spiral  configuration  is  still  an 
imanswered  question.  Experiments  and  further  modeling  to  address  the  high  velocity 

dependence  of  //(v)  are  slated  to  take  place  in  the  near  future. 
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